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Abstract 

In this article, we present analytical expansion results of two single mass scale four-loop vacuum 
integrals in d = 3 — 2e dimensions. After finding hypergeometric representations with half- 
integer coefficients, we use algorithms which we implemented in FORM to expand these in 
terms of nested sums. 
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1. Introduction 

The series of the pressure of QCD in weak coupling expansion is not analytical in the 
coupling (7^. The reason for this is that hot QCD is a multiscale system. One way of 
dealing with this problem is using effective field theory approach [l|2|3j . In this approach 
the different scales are separated within different effective theories. At temperatures well 
above Tc all the relevant information of QCD is contained within a three dimensional 
effective field theory, the so-called EQCD. For the specific problem mentioned above, 
EQCD has the following set of master integrals up to four loops [4], obtained using 
methods from |5|6|7j : 




(1) 



The full/dotted lines stand for a massive/massless scalar propagators. These Masters 
need to be computed. There are several methods for doing this; among the most frequently 
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used are solving the corresponding difference equations [7], as well as the Mellin-Barnes 
method |8l9j . Here we will only be using Mellin-Barnes method for obtaining hypergeo- 
metric representation of some Masters, while as a second step we use nested sums |10lllj 
for expanding these hypergeometric functions analytically around d = 3 — 2e dimensions. 



2. Mellin-Barnes method 



The Mellin-Barnes method for evaluating massive Feynman integrals consists of writing 
down a massive propagator in terms of a massless one [5] : 
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where the contour of integration is chosen such that the poles of r(— s) arc to the right 
of the contour. One simple example is the one-loop integral of eq.(3), which using the 
Mellin-Barnes method gives us a massless one-loop integral which in turn can be evaluated 
in terms of Gamma functions. At the end one can perform the complex Mellin-Barnes 
integral closing the contour to the right, obtaining: 
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What is left to do now is the sum over Gamma functions. Fortunately, there are 
methods developed [11112] for expanding in e such hypergeometric sums in terms of 
objects called nested sums. Note that in the specific example above, for d = 4 — 2e one 
would have integer valued coefficients in F-functions, while for d = 3 — 2e one is left with 
half-integer valued coefficients. 

3. Nested Sums 

S-Sums [lOlllj are defined recursively as: S{n) = 1 for n > , S'(n) = for n < and 

n ^ 

S{n;mi, . . . ,mfc;xi, . . . ,a;fc) = ^ T^S'(i;m2, . . . ,mk;x2, .. .,Xk) ■ (4) 
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Gamma functions can be naturally expanded in terms of S-sums, according to: 

^^(i^ = r(., + i).xp(-|:.'H)!50;M)). (5) 

So if we can represent Feynman integrals in terms of one infinite sum over Gamma 
functions, we can expand the Gamma functions in terms of S-Sums, use their algebra 
[lOlllj to reduce products of S-sums to a linear combination of S-sums and finally do 
the last infinite sum over S-sums. This last sum gives S-sums with argument infinity. 
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which are, per definition, muhiple polylogarithms (MPL)[T3j. In eq.(6,7) MPL simphfy 
further to multiple C- values [14 . There are packages for the expansion of hypergeometric 
functions with integer valued coefficients [15116117] and of some classes of functions with 
half-integer coefficients [18]. However, for some of our Masters, in particular the one of 
eq.(6), these are not sufficient, what prompted us to implement the relevant algorithms 



4. Applications 

Here we apply the methods discussed in previous sections to two four-loop single mass 
scale Feynman integrals from EQCD. We first apply the Mellin- Barnes method and get 
hypergeometric representations. These consists in the case of eq.(7) of one hypergeomet- 
ric function, which in d = 3 — 2e is given by 4F3 1^1, 1 — e, 2e, i -I- e; 2 — 2e, 1 -I- e, | — e; 1^ 
and in case of eq.(6) of nine different hypergeometric functions, of which the highest one 
is (^1, 1 - e, I -f e, I -t- e, 2 + 2e; 2 - 2e, I - e, I + e, 2 + e; 1^ . As mentioned above, we 
implemented the algorithms for expansion in e of half-integer valued balanced hyperge- 
ometric functions in FORM [19] . Hypergeometric representations in d dimensions and 
expansions in d = 4 — 2e will be published elsewhere [20] . Note that both integrals are 
normalized by the appropriate power of the massive one-loop tadpole integral. 



of [12]. 
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where = £ E = 0.987441... and a, = Li,(i). 



5. Summary 

We implemented in FORM algorithms for expanding hypergeometric functions, in- 
cluding expansion around half-integer balanced values. As an application we showed 
expansions in c? = 3 — 2e of two four-loop bubbles needed in the context of hot QCD. 



References 

[I] T. Appelquist and R. D. Pisarski, Phys. Rev. D 23 (1981) 2305. 
[2] P. H. Ginsparg, Nucl. Phys. B 170 (1980) 388. 

[3] E. Braaten and A. Nieto, Phys. Rev. D 51 (1995) 6990 

[4] Y. Schroder, Nucl. Phys. Proc. Suppl. 116 (2003) 402 

[5] K. G. Chetyrkin and F. V. Tkachov, Nucl. Phys. B 192 (1981) 159. 

[6] F. V. Tkachov, Phys. Lett. B 100 (1981) 65. 

[7] S. Laporta, Int. J. Mod. Phys. A 15 (2000) 5087 

[8] E. E. Boos, A. I. Davydychev, Theor. Math. Phys. 89 (1991) 1052 [Teor. Mat. Fiz. 89 (1991) 56]. 
[9] V. A. Smirnov, Springer Tracts Mod. Phys. 211 (2004) 1. 
[10] J. A. M. Vermaseren, Int. J. Mod. Phys. A 14, 2037 (1999) 

[II] S. Moch, P. Uwer and S. Weinzierl, J. Math. Phys. 43 (2002) 3363 
[12] S. Weinzierl, J. Math. Phys. 45 (2004) 2656 

[13] A. B. Goncharov, Math. Res. Lett. 5 497 (1998) 

[14] J. M. Borwein et al., Trans. Am. Math. Soc. 353 (2001) 907 

[15] S. Moch and P. Uwer, Comput. Phys. Commun. 174 (2006) 759 

[16] T. Huber and D. Maitre, Comput. Phys. Commun. 175 (2006) 122 

[17] S. Weinzierl, Comput. Phys. Commun. 145 (2002) 357 

[18] T. Huber and D. Maitre, Comput. Phys. Commun. 178 (2008) 755 

[19] J. A. M. Vermaseren, arXiv:math-ph/0010025| 



4 



[20] E. Bejdakic and Y. Schroder, in preparation 



5 



